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SIZE BIAS AND DIFFERENTIAL LENSING OF STRONGLY LENSED, DUSTY GALAXIES IDENTIFIED IN 

WIDE-FIELD SURVEYS 
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ABSTRACT 

We address two selection effects that operate on samples of gravitationally lensed dusty galaxies 
identified in millimeter- and submillimeter- wavelength surveys. First, we point out the existence 
of a "size bias" in such samples: due to finite source effects, sources with higher observed fiuxes 
are increasingly biased towards more compact objects. Second, we examine the effect of differential 
lensing in individual lens systems by modeling each source as a compact core embedded in an extended 
diffuse halo. Considering the ratio of magnifications in these two components, we find that at high 
overall magnifications the compact component is amplified by a much larger factor than the diffuse 
component, but at intermediate magnifications (~10) the probability of a larger magnification for the 
extended region is higher. Lens models determined from multi- frequency resolved imaging data are 
crucial to correct for this effect. 

Subject headings: gravitational lensing — galaxies: luminosity function, mass function — galaxies: 
abundances — methods: numerical 



1. INTRODUCTION 

Surveys of the extragalactic sky at millimeter and sub- 
millimeter wavelengths have identified an optically dim 
galaxy population in which rapid, obscured star forma- 
tion powers far infrared (FIR) luminositie s occasionally 
exceeding IQ-^^ L^:^ (e.g., ISmail et aTl 19971: iBarger et all 
[19981 : iHughes et all [19981 : lEales et al.l[T999[ ). Surveys of 
these sources at 850 um using SCUBA were confusion- 
limited, with several thousand sources per square degree 
at mJy- level fiuxes (Coppin et al. 2006). The source 
density falls steeply with increasing fiux: at 850 /im 
there are only a h andful of sources p er square de- 
gree above 10 mJy (jCoppin et al.l [20061 ). In the first 
millimeter- wavelength survey to cover tens of squ are de- 
grees with similar sensitivity, IVieira et al.l ()2010l ) found 
that the rarest objects (< ldeg~^) are far more abun- 
dant than would be expected from the extrapolation of 
lower fiux density sources. Such a population is a nat- 
ural consequence of the gravitational lensing of intrinsi- 
cally fain ter galaxies by interv ening galaxies and galaxy 
clusters ( Negrello et al. "2007'; "Hezave h fc Holder! 120111 
iBethermin et al. 2011), and thus the selection of very 
bright submillimeter sources from large area surveys has 
proven t o be an efficient way to identify highly-magnified 
objects (jNegrello et al.ll2010[ ). 

The gravitational magnification of these sources al- 
lows them to be examined in more detail: they are 
seen at improved source-plane resolution (due to the 
magnification from gravitational lensing), allowing the 
study of the star for mation process on the s cale of gi- 
ant molecular clouds (jSwinbank et al.l l2010 a^ ; further- 
more, the overall increase in fiux allows higher signal- 
to-noise measurements of the spectra l energy distribu- 
tion (SEP) and molecular lines (Riecherset al.l 
Lestrade et al.l 120111 : iFraver et al.ll2011l : iCox et^aP 
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the gas conditions. 

Because gravitational lensing is a geometric effect, a 
given position in the lensed object will be magnified 
equally at all wavelengths. However, the finite extent 
of the background galaxy will result in variations in 
the magnificati on applied to differen t regions within the 
galaxy (Blandf ord fc NaravanI ri992[ ) : the observed SED 
(e.g., Blain 1999), as well as ratios of spectral lines (e.g., 
jDowne s et al. 1995), will be distorted by this differen- 
tial magnification if there are spatial variations in the 
physical conditions within the source. 

The interpretation of observations of gravitationally 
magnified submillimeter galaxies (SMGs) are necessarily 
comp licated by this possibility of differential magnifica- 
tion. iBlainI ([l99 9) explained the excess mid-IR emission 
in some lensed SMGs through preferential magnification 
of compact hot regions, which can fiatten the observed 
spectrum at wavelengths shorter than the SED peak and 
thereby increase the inferred temperature and luminos- 
ity of the source. As shown below, a variety of effects 
are possible, with the possibility of preferentially magni- 
fying either compact or diffuse regions, depending on the 
relative position of source and lens. 

The sensitivity advantages associated with observing 
lensed sources will make these galaxies prime targets for 
studying the physics of SMGs. However, biases intro- 
duced by the selection of lensed galaxies must be con- 
sidered carefully when extrapolating physical properties 
(e.g., temperatures, luminosities, sizes) determined from 
these samples to the SMG population as a whole. 

In this paper we consider two effects that operate 
in lensed SMG populations selected by their high fiux 
in millimeter and submillimeter surveys. First, we use 
population models for S MGs, previously published in 
iHezaveh fc Holderl ()2Q11[ ). to examine the total magni- 
fication present in such samples and point out the possi- 
bility of a strong size bias that such selection introduces. 
Second, we examine the differential magnification of a 
two-component galaxy model as a function of properties 
of the source, lens, and the source- lens alignment. 

For all calculations below, we assume a spatially fiat 
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universe with Qrn = 0.222, h = 0.71, n^ = 0.96, and 
(78 = 0.801. 

2. MODEL PARAMETERS 

2.1. Population Modeling 

iHezaveh fc Hold er (2011") reproduced the observed 
abundance of bright SMGs at 1.4 mm b y combining 
models of the redsh ift distribution of SMGs (jLacev et al.l 
I2Q1QI : iMarsden et al. 2011) and a pop ulation of ellipti- 
cal len ses with masses described by the iSheth fc TormenI 
()1999f ) mass function. For our exami nation of the popu- 
l ation of lensed SMGs, we adopt the IHezaveh fc Holderl 
([2011) model, with the underlying S MG popul ation de- 
scribed by the semi-analytic model of ( Lacev et a l. 2010). 
The lensing magnification probability distribution is cal- 
culated as a function of source size by numerically cal- 
culat ing the lensing cros s -secti on for a particular source 
size. 'Hezav eh fc Holderl (J2Q11I ) assumed a mean source 
size and hence used a single lensing probability which 
was applied to the number counts. In this work we are 
examining the effect of the natural size distribution of the 
SMGs on the observed number counts and hence apply 
four lensing probability functions corresponding to four 
different source sizes. We simply assume a source popu- 
lation consisting equally of point sources (small enough 
to avoid magnification damping effects at the 1-100 range 
for galaxy mass lenses), and sources with radii of 1.0, 3.0, 
and 8.0 kpc and study their apparent relative abundance 
after lensing. Such sizes have b een measured in resolved 
SMGs re.g.. lTacconi et al.l[2008l) . 

2.2. Lens Models 

For the purpose of studying the effects of differential 
lensing of a single lens system we use our ray-tracing 
code to simulate lensed images of extended sources. We 
model our lens as a Singular Isothermal Ellipsoid defined 
by its mass, ellipticity and angle and place the lens at 
Zd = 0.5. The source is modeled with a slightly ellip- 
tical morphology as a two component model consisting 
of an extended component and an inner core with the 
ratio of the outer radii defined as aR = Rcore/Rext- The 
unlensed fiux ratio of the two components is defined as 
Pf = Fcore/Fext = 1 and the total fiux Ft is simply 
the sum of the two component fiuxes. After placing this 
source behind a lens the total magnification /ix is defined 
as the ratio of the lensed to unlensed fiux of the source 
and /3p is the observed flux ratio of the two components. 
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This morphology is motivated by the star forming 
structures observed in galaxies, which typically include 
large reservoirs of cold dust (extended component) in ad- 
dition to small clumps o f hot dust (core) in the vicinity 
of sta r forming regions (jPohlen et al.l 120101 : iHaan et all 
1201 If ). Although real galaxies may host multiple clumps 
with spatial offsets relative to the cold dust, this simple 
model can demonstrate the systematic effects caused by 
differential lensing of more complex sources. 
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Fig. 1. — Lensed number counts of 1. 4 mm (220 GHz) sources 
based on Durham semi-analytic m odel (| Lacev et al. 2010, as de- 
scribed in 'Hezaveh &: Holdeiil20l T) assuming four different source 
radii. The solid lines are calculated for a lens population with an 
ellipticity (1-axis ratio) of e = 0.1, while the dot-dashed lines are 
for e = 0.3. The lensing mass distribution is renormalized for the 
two cases so that the SPT (IRAS-excluded) dusty galaxy source 
counts (black circles; Vieira et al. 2010) are matched. The bottom 
panel shows the ratio of the lensed counts of each population to the 
total lensed number counts showing that at higher observed fluxes 
there is a strong bias towards selecting the most compact sources 
from the input population. 




Fig. 2. — From the source count model of Figure ^ the distri- 
bution of source magniflcations in the population of SMGs at two 
different observed flux densities. Green, blue, red, and cyan curves 
(bottom to top) correspond to point sources, and sources with radii 
of 1, 3, and 8 kpc, respectively. Solid curves show the magniflca- 
tion distributions at observed flux cut of 25 mJy and the dashed 
curves show the distribution at 2.5 mJy. For typical SMG spec- 
tra, these observed flux densities roughly correspond to 350 and 
35 mJy for Herschel surveys selecting sources at A=500 /xm. Un- 
magnifled sources are excluded. The high magniflcation sources are 
an increasing fraction of the sources at the highest flux densities. 
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Fig. 3. — Cumulative magnification cross-sections of different source morphologies (insets) for a 7 x 10^^ Mq galaxy (mass within the 
Einstein radius) at z = 0.5. The source is placed at Zs = 2.0. The sizes in the legend correspond to the radius of the star forming clumps. 
All the insets are on the same scale. Left: Source models with varying numbers of clumps, to fill the same sky area as the 1 kpc source. 
Right: Source models with a fixed number of clumps in fixed locations. 
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3. RESULTS 

We applied the models to investigate two related but 
separate effects: 1) given some dispersion in intrinsic 
source sizes, how does selecting the brightest lenses af- 
fect the inferred source size distribution; 2) for a sin- 
gle multi-component source (modeled here as one small 
emission region embedded in a larger emission region), 
how does lensing selection affect the relative amount of 
flux received from regions of different sizes? 

3.1. Size Selection 

The modeling shown in Figure [1] demonstrates that the 
distribution of source sizes derived from samples of lensed 
objects are unlikely to be representative of the underly- 
ing source size distribution. Starting from a population 
of sources that are equally distributed among four sizes 
(Section 12. ip , we find that the selection of millimeter- 
bright sources (> 10 mJy at 1.4 mm) introduces a strong 
selection bias toward the most compact sources. Com- 
pact sources can therefore be expected to be overrepre- 
sented by a factor of two or more in a survey with a 
high flux density threshold, while larger sources would 
be nearly completely missed. 

This result is not trivial, as it represents the balance 
between the higher total magnification achievable for a 
compact source and the lower probability of aligning the 
compact source with a lensing caustic to achieve high 
magnification. Figure [1] demonstrates that the bias per- 
sists across a range of realistic lens ellipticities, with a 
stronger bias toward compact objects for more elliptical 
lensing potentials. The ellipticities of real lenses will be 
distributed through and beyond the e = 0.1 — 0.3 range 
shown, and this unknown distribution function would be 
required in order to predict a robust correction for the 
size bias in observed populations. 

The preference for compact sources at the highest flux 
densities can be explained by examining the distribution 
of magnification in the lensed population as a function 
of observed flux density (Figure Ej). At larger observed 
flux densities, the distribution of source magnifications in 
the lensed sample (unlensed sources are excluded) shifts 
to larger values. The distribution of magnifications is 
dependent on many aspects of the source population 
model and the distribution of lensing halo shapes and 
masses. However, the basic conclusion that sources of 
the highest observed flux density are dominated by high- 
magnification sources is robust. The brightest sources 
are drawn from a wide range of intrinsic fluxes or, equiv- 
alent ly, objects of the same observed flux will derive from 
a wide range in /i. A sample of sources with a high flux 
density threshold will best identify very high magnifica- 
tion sources, which also provides the highest source-plane 
spatial resolution for studying the ISM on small scales in 
these distant galaxies. 

We also study the possible effect of the dumpiness of 
the source light profile on lensed number counts. High- 
spatial resolution observations have shown that lumi- 
nous SMGs are compo sed of a few star forming knots 
(|Swinbank et al.ll2010b[ ). Since each knot is smaller in 
size than the overall extent of the galaxy, it is possible 
that the probability of achieving high magnification will 
change significantly with the size of the knot. To study 
this effect we calculate the magnification cross-section. 




Fig. 4. — A lens model showing the effects of differential mag- 
nification. The total magnification of the three subcomponents 
are 11.7, 17.1, and 5.2 for the extended source (gray) compact 
source inside the caustic (blue), and the doubly imaged source 
(light green) respectively. Here the second compact component is 
assumed in order to emphasize the strong dependence of magnifi- 
cation on the exact spatial configuration of the components. 

the solid angle in the lens plane over which the magnifi- 
cation will exceed a given value, for a source of roughly 
fixed size (1 kpc radius) but broken into clumps of five 
different sizes, radii of 1 kpc, and 500, 200, 100, and 50 
pc. 

In Figure [3l we compare the lensing cross-section for 
the various source morphologies. In the left panel, where 
for each clump size we distribute them over the entire 
region covered by the 1 kpc source, we find that the 
dumpiness of the emission has very little impact on the 
overall lensing magnification. The clumps sample the 
same region of the lens as the single source, so the mag- 
nification averaged over the ensemble of clumps is similar 
to that for the single source. In the right panel, we di- 
vide the 1 kpc source into just six clumps of different 
sizes. The sparser sampling of the lens leads to small 
but measurable disparities in the magnification cross sec- 
tion between the models with different clump sizes. The 
smallest clumps have the largest possible magnifications, 
but also have a reduced cross section between /i = 10 
and 20 for this combination of clump and lens geom- 
etry. Together, these tests demonstrate that the total 
extent of the background galaxy, rather than its division 
into clumps, is primarily responsible for determining how 
likely it is to be lensed, with the small scale structure 
having some impact on the total magnification. 

3.2. Differential Magnification 

The selection of high-magnification sources assures 
close alignment between the lensing caustics and regions 
of bright emission on the SMG. This increases the im- 
portance of the differential magnification effect, which 
is demonstrated in Figure HI Here, adjacent compact re- 
gions that straddle the tangential caustic differ in magni- 
fication by a factor of three, while the total magnification 
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Fig. 5. — Top left: Magnification ratio as a function of position in the source plane, with e = 0.3, q;r = Rcore/Rext = 1/8, and intrinsic 
flux ratio (3-^ = 1] this is our fiducial model. The total magnification is indicated in all panels with white contours. The red countour traces 
the line of constant total magnification that passes through the minima in the magnification ratio (/5p). The source model is outlined in 
the lower right. Top right: The fiducial model, but with lens ellipticity decreased (e = 0.1). Bottom left: The fiducial model, but with the 
source geometry altered (q;r = 1/4). Bottom right: The fiducial model, but with the source size doubled. 
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of the extended component that surrounds them is com- 
parable to the magnification of the brightest component. 
For real sources that are a composite of multiple com- 
ponents of significantly varying conditions, the effects of 
differential magnification on the observed SED and spec- 
tral line ratios can be dramatic (Serjeant 2012) . 




Fig. 6. — Normalized distribution of flux ratios at each magniflca- 
tion. At intermediate magniflcations (e.g. ji ~ 10) there is a higher 
probabihty of a larger magniflcation for more extended objects. At 
higher magniflcations the situation is reversed. This combined with 
the overall magniflcation distribution at an observed flux (Figure 
[2)) results in a serious difliculty in interpreting the intrinsic flux 
ratio of different components. 

Starting from the simple source structure described in 
Section 12. 2[ we consider lensing effects for several galaxy 
and lens variations in Figure [51 A fiducial model with 
e = 0.3, aR = 1/8, and intrinsic fiux ratio /3f = 1 is 
shown in the top left. The source plane is shown, the 
color scale represents the (logio) observed fiux ratio, cor- 
responding to the ratio of magnifications (/3p) of the two 
components as a function of the source- lens alignment. 
The ratio varies by a factor > 10 across this image, in- 
cluding substantial area over which the extended compo- 
nent can be accentuated by a factor of several. 

The lens and source parameters both strongly affect 
the observed ratios. The top left and right panels show 
that in a less elliptical lens, the magnification ratio varies 
in a smaller region. The similarity between the size of 
the diamond caustic and the "core" emission component 
results in less variation in the magnification ratio across 
the source. The left panels demonstrate that making the 
sizes of the two components more similar decreases the 
variations in magnification ratio. Increasing the size of 
the whole source (lower right) slows the rate of change 
of the magnification ratio. The red contour in each plot 
shows the curve of constant /ix that connects the points 
of minimum /3p. This curve is always very close to the 
caustic, where slight changes in source- lens alignment 
can mean dramatic changes in magnification ratio. 

The probability distribution for /3p at fixed /ix is shown 
in Figure [6l for the model in the top left panel of Figure [51 
At moderate magnifications, the extended component is 
often magnified by a larger factor, while the highest mag- 
nifications over emphasize the core but with a large range 
in amplitude. Considering Figures [2l and [5l together, we 



see that for a sample selected with a high fiux density 
limit, the distortion of the SED cannot be assumed to be 
understood. 

A more optimistic view of the situation can be taken 
from Figure [5l A good lens model can greatly reduce 
the uncertainty in the magnification distribution across 
the source. Regions where the fiux is dominated by the 
diffuse component are confined to those areas just out- 
side the "diamond caustic" that are well-separated from 
the cusps. High-resolution imaging with ALMA or Hub- 
ble Space Telescope will provide lens models that will be 
useful for these purposes. While it would be preferable to 
have high-resolution imaging at all frequencies, for sim- 
ple questions about whether the fiux is dominated by the 
compact or diffuse regions it could be sufficient to have 
a good lens model at a single frequency. More advanced 
inferences are available with multiple wavelengths of re- 
solved imaging or by assumin g a physically motivated 
size structure (e.g.. lBlainlll999[ ). 

4. CONCLUSION 

Strongly lensed SMGs permit close examination of 
these rapidly star-forming galaxies with short integra- 
tion times. However, the selection of lensed objects and 
the strong lensing itself introduce biases that must be 
considered when interpreting the observed galaxy prop- 
erties. 

We have demonstrated that strongly lensed SMGs se- 
lected by a flux cut in wide surveys are affected by a "size 
bias". They are more likely to include more compact 
galaxies, which correspond to higher surface brightness 
sources at a given unlensed flux. More speciflcally, as 
one probes samples with higher observed flux the distri- 
bution of the size of the sources is increasingly skewed 
toward more compact objects. 

We have also studied the effects of differential lensing 
of SMGs. We model the SMGs as having two compo- 
nents, a compact core embedded in a more extended re- 
gion, and measure the relative flux of these components 
in the presence of strong lensing. We flnd that the rela- 
tive gain for compact and extended components depends 
sensitively on the source-lens alignment. At high total 
magniflcations the observed emission from the compact 
core is amplifled by a larger factor relative to the emis- 
sion from the diffuse region. At intermediate magniflca- 
tions (e.g., /i ~ 10), the situation is often reversed with 
the diffuse emission being magnifled by a larger factor. 
A similar range of effects can be expected for spectral 
lines. 

This suggests that without applying the required cor- 
rections from lens models, the SED and the ratios of 
molecular lines could be signiflcantly biased. Fortu- 
nately, the relative magniflcation of the compact and ex- 
tended regions depends on the geometry of the caustics 
and the source-lens alignment in a highly regular man- 
ner. This effect can therefore be well understood by con- 
structing reliable lens models for each source, which will 
require spatially resolved imaging data. 
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